Electron-phonon effects and transport in carbon nanotubes 
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We calculate the electron-phonon scattering and binding in semiconducting carbon nanotubes, 
within a tight binding model. The mobility is derived using a multi-band Boltzmann treatment. At 
high fields, the dominant scattering is inter-band scattering by LO phonons corresponding to the 
corners K of the graphene Brillouin zone. The drift velocity saturates at approximately half the 
graphene Fermi velocity. The calculated mobility as a function of temperature, electric field, and 
nanotube chirality are well reproduced by a simple interpolation formula. Polaronic binding give 
a band-gap renormalization of ~70 meV, an order of magnitude larger than expected. Coherence 
lengths can be quite long but are strongly energy dependent. 

PACS numbers: 73.63.Fg 
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Carbon nanotubes have enabled novel electronic de- 
vices, including quasi-one-dimension field-effect transis- 
tors and electro-optical devices 0, ll| • Much effort has 
gone into determining the transport properties, which are 
crucial for these and other applications. Nanotubes ex- 
hibit high mobility at low electric fields, exceeding the 
best semiconductors at room temperature 0. At high 
fields, the mobility is dramatically reduced by optical 
phonon emission, leading to interesting saturation be- 
havior 0,0,0- However, little is known about how the 
mobility varies with temperature, electric field, and tube 
diameter and chirality. 

Here we calculate the electron-phonon interactions and 
drift velocity in semiconducting nanotubes, within a 
standard tight-binding approach. Our results are con- 
sistent with available measurements, and provide a de- 
tailed microscopic picture of phonon scattering. For ex- 
ample, we find that there is a broad maximum in coher- 
ence length at intermediate energies, centered at 100-150 
meV, with coherence lengths of over a micron at room 
temperature, and much larger at low temperature. At 
large fields, the mobility is limited primarily by inter- 
band scattering from the LO phonons corresponding the 
corners of the graphene Brillouin zone. 

At the same time, we provide a broad overview of the 
transport, by calculating the mobility over a wide range 
of temperature, electric field, and nanotube structure. 
The low-field mobility depends strongly on nanotube di- 
ameter and temperature, while at high field the drift ve- 
locity approaches a roughly universal value. We show 
that these dependences can be well described by a sim- 
ple formula. 

In addition, we find that polaronic binding gives a sig- 
nificant bandgap renormalization, reducing the gap by 
around 70 meV over a range of tube diameters. This is 
more than an order of mEignitude larger than suggested 
by previous calculations [a 1^, and is particularly im- 
portant for larger-diameter tubes, where it represents a 
significant fraction of the bandgap. 

Our calculations use a standard tight-binding descrip- 



tion for the electronic structure (IQl] . The phonon en- 
ergetics are described using an improved model some- 
what similar to that of Aizawa et al. ■ We model the 
electron-phonon (e-ph) interaction by the Su-Schriefer- 
Heeger (SSH) model [12|, with matrix element t — to ~ 
gSu dependent on the change of the nearest neighbor C- 
C distance (Su), where to — 5 eV. We take the e-ph 
coupfing constant to be g = 5.3 eV/A as pre dicted theo- 
retically for a related molecular problem [13j, consistent 
with fits to the Peierls gap in conjugated polymers [l^. 
The Fourier transformed SSH Hamiltonian here is 
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where M^^^ c>c gN is the e-ph coupling IF]; u\^^ 
(vk) denotes creation (annihilation) of an electron in the 
conduction (valence) band; a^^^ is a phonon creation 
operator with wavevector ~q and phonon band index fi — 
1...6; N is the number of the primitive unit cells each 
containing two carbons. The indices k and q each label 
both the continuous wavevector along the tube axis and 
the discrete circumferential wavevector. 

Since conduction and valence bands are symmetric in 
carbon nanotubes, we consider a filled valence band with 
a single electron in the conduction band. The behavior 
of a single hole would be identical. The charge carrier 
with wavevector k and energy Sk can be scattered to a 
state with wavevector k + q by absorbing (emitting) a 
phonon with wavevector q {—q), and having energy fkUq 
(= huj^q) such that the net momentum and energy are 
conserved. 

Scattering and binding correspond to the imaginary 
and real parts of the self-energy. We calculate these in 
the Random Phase Approximation. The scattering rate 
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for an electron of wavevector k is 
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where is the scattering time, and the phonon occupa- 
tion number Uq is given by the Bose-Einstein distribution. 
The energy shift due to electron-phonon couphng is 
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in the Hmit i5 — > 0. 

The results for the binding energy and the scattering 
rate of a charge carrier in the first band are shown in 
Fig. ^for a tubes of diameter d=2.0 nm. The binding 
energy in Fig. ^ is nearly independent of temperature 
over the entire energy range, and only weakly energy- 
dependent at low energy. The binding energy increases 
sharply at resonance with the optical phonon K^o- The 
scattering rate (inverse lifetime) in Fig.^ shows a strong 
temperature dependence at low energy, but is two orders 
of magnitude larger and nearly independent of tempera- 
ture at high energy. 

Virtually all of the scattering is by phonons corre- 
sponding to a small region near the F and K points of the 
2D graphene Brillouin zone, and Fig. ^ shows the con- 
tributions of the respective bands to the scattering. At 
low electron energy, the scattering is by acoustic phonons 
near the F-point. The lowest-energy phonon band gives 
negligible scattering and is not shown. The second band 
is the acoustic phonon contribution which peaks at the 
van- Hove singularity (bottom of the conduction band). 
The third band is the radial breathing mode which gives 
a peak labelled Frbm at 15 meV in (25,0) tube. The 
next phonon mode with a non-negligible e-ph coupling is 
the longitudinal acoustic phonon at K-point which gives 
a weak scattering peak labelled Kla- 

By far the strongest coupling is to the LO phonons, 
which give scattering roughly two orders of magnitude 
stronger than the acoustic modes. The Klo mode corre- 
sponds to inter-band scattering. It is the most important 
mode at high field, both because it is strongest, and be- 
cause it can scatter electrons before they reach the energy 
of the Flo mode. 

In Fig. nji we show the coherence length, defined as 
Lk = VkTk, where Vk = dsk/dk is the band velocity. 
(We approximate this by the bare velocity Vk — dsk/dk, 
which we calculate to be the same to within a few per- 
cent except near resonance with the optical phonons Klo 
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FIG. 1: (color online). Phonon scattering and binding vs. 
electron energy (relative to conduction band edge £c), for 
a (25,0) tube (d=2.0 nm): (a) binding energy, Eq. @; (b-c) 
inverse lifetime, Eq. Q; and (d) coherence length. The en- 
ergy range shown includes only the first band. In (a) , (b) , and 
(d), solid red curves are for T=300 K, dotted blue curves are 
for T=10 K. Note the large change in scale between low and 
high energy, x50 in (b) and (c), and x500 in (d). (c) Decom- 
position of e-ph scattering rate for 300K: total (red curve); 
F-point contributions (green curves); K-point contributions 
(blue curves). r_H_Bi\/ is radial breathing mode. The inset 
shows 2D graphene Brillouin zone, with the integration area 
shown by the circles of radius equal to 0.1 of the F-M distance, 
which contributes virtually 100% of the total scattering; lines 
are only schematic. 



and Flo-) Note that the coherence length is strongly 
energy-dependent, with a broad maximum between the 
breathing mode and the optical phonons. It is also quite 
sensitive to temperature, but even at room temperature 
the coherence length can be well over a micron, consistent 
with lengths inferred from experiment 0,0- Charge car- 
riers can be injected into nanotubes well above the band 
edge T^, suggesting the possibility of ballistic or even 
quantum-coherent devices over a range of length scales. 
At higher energy, the carriers have a temperature inde- 
pendent coherence length of around 20-40 nm, depending 
on energy. This is consistent with the lengths inferred 
from experiments for metallic tubes [a, |a, [Tf . 

For applications such as device modeling, what is 
needed is the total effect of scattering on the transport. 
This is expressed as a mobility, which we calculate by 
solving the steady-state multi-band Boltzmann equation 
in the presence of an electric field. The electron momen- 
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FIG. 2: (color online), (a) Drift velocity vs. electric field 
at T=300 K (red curve) and T=10 K (black curve) in (25,0) 
tube, (b) Inverse mobility at T=300 K vs. i?-field; straight 
lines show corresponding fits to Eq. with Vs— 5.0x10^ . 
Fitted values of /lo are /io=65000 cm'^/Vs for (25,0) tube (red 
triangles); /io=35500 cm^/Vs for (19,0) tube (green squares); 
and ^0 = 15000 cm^/Vs for (13,0) tube (blue circles). 

turn distribution function is 

eEdgk ^-^ . , . 
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where Wk,k+q is given in Eq. lO. Using the non- 
equihbrium distribution function gk from Eq. |0J, we cal- 
culate the drift velocity v and the mobility [i: 

- = E5'c^=Mi? (5) 

k 

The drift velocity versus i?-field is shown in Fig.|2jL, for 
a (25,0) tube. We show results at room temperature and 
at low temperatures, in the limit of low carrier density. 
The drift velocity saturates for fields E > 0.5 V//xm. At 
low temperature we find negative differential resistance, 
i.e. drift velocity decreasing with increasing field. A sim- 
ilar observation was reported previously [l7|, and it was 
speculated to result from population of the second band 
at high field. However, we find that the behavior persists 
even if we restrict the basis set to include only the first 
band. 

The inverse mobility shown in Fig. can be fitted 
rather well by a simple linear function of the i?-field, 
over a wide range of applied E-ficlds: 

« Mo 1 + v;'E, (6) 
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FIG. 3; (color online), (a) Zero-field mobility vs. tube di- 
ameter, for tubes of many different chiralities. Temperatures 
are T=450 K (red), T=300 K (green), and T=150 K (blue). 
Values are based on linear extrapolation of /i^^, Eq. 0, for 
i?-field between 0.01 and 0.5 V / fim. Black curves are fitted to 
Eq. Q. (b) Polaron binding energy (black squares) vs. tube 
diameter. Red circles show portion of binding energy from 
optical phonons, and blue triangles show the much smaller 
binding energy (xl5) from acoustic phonons. 

where is the zero-field mobility and is the satura- 
tion velocity. This is analogous to the linear dependence 
R = Rq + V/Iq found for the dependence of the resis- 
tance R of metallic nanotubes with applied voltage V 
1^], where Iq is a saturation current. The saturation drift 
velocity Vs « 5.0±0.3 x 10^ cm/s is nearly independent of 
the tube diameter and temperature. The lack of strong 
dependence on diameter is perhaps surprising, but this 
velocity is intriguingly close to half of the maximum band 
velocity (the Fermi velocity of graphcne, 9.8 x 10'' cm/s 
for the value of the tight-binding matrix element t used 
here). 

One of the greatest obstacles to fabricating repro- 
ducible nanotube devices, is the difficulty of controlling 
(or even knowing) the tube structure. This structure, the 
diameter and chirality, is specified by two indices (m,n). 
In Fig. we show how the zero-field mobility depends 
on tube structure. We find that the dependence on diam- 
eter and temperature can be described by a single simple 
empirical relation: 

300K fay 

1 \ Inm / 

with fii = 12000 cmVVs and a = 2.26. (The arbitrary 
constants 300K and 1 nm are simply to give fii units of 
mobility, independent of a.) The structure dependence is 
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FIG. 4: (color online). Steady-state distribution function for 
an electron in a (25,0) tube at T=300 K. E-fields (from bot- 
tom to top on right side) are E=0.001 (black), 0.01 (blue), 
0.03 (light blue), 0.1 (green), and 0.5 (red broken curve) 
W / fj,m. Distributions are shown with respect to (a) energy, 
and (b) wavevector along tube axis. (Curves include Gaussian 
broadening of 2 meV.) The top axis in (b) shows correspond- 
ing energy. 



not entirely captured by the diameter, nor is the 1 /T tem- 
perature dependence very precise, but overall the agree- 
ment is quite good. Thus the dependence of mobility on 
temperature, field strength, and even nanotube structure 
can be captured to useful accuracy (~10%) vifith a simple 
3-parameter description, Eqs. (|6I7|I . 

It might seem surprising that, while the scattering in 
Fig.^is strongly energy-dependent, the mobility in Fig. [3 
varies smoothly with increasing field. This can be un- 
derstood by considering the distribution function [or 
equivalently f/s] obtained by solving the Boltzmann equa- 
tion. This is shown in Fig.01 At very low field, the scat- 
tering is entirely by acoustic phonons, giving a smooth 
distribution in k. (The distribution in energy refiects the 
van Hove singularity in the density of states.) With in- 
creasing field, the electron has an increasing probability 
of reaching the optical phonon energy, at which point it is 
quickly scattered. This leads to a step in the distribution 
near the phonon energy (0.18 eV). With increasing field, 
the distribution becomes flatter below this energy, until 
at high field the distribution is almost fiat in k, with a 
sharp step at the phonon energy. From Eq. ((^Jl, we might 
expect the crossover from acoustic to optical phonons as 
the dominant scattering mechanism to occur roughly at 
fields E ~ Vs/no, corresponding to '^0.08 V//zm in Fig.^ 
This is indeed roughly the field strength where interme- 
diate behavior is seen in gk- Then from Eqs. (|5I7|I . we 
can anticipate that the crossover field varies with tem- 



perature and tube size approximately as d "T. 

Finally we return to the polaronic energy shift. Fig- 
ure shows the shift calculated with Eq. as a func- 
tion of nanotube diameter. (This refers to the band-edge 
states, Ek—Sc in Fig. ^.) We find binding energies of 
around 35 meV, with only a weak dependence on diam- 
eter over the typical range. This corresponds to a band- 
gap renormalization of 70 meV due to electron-phonon 
interactions, which is quite significant for large-diameter 
tubes. This is more than an order of magnitude larger 
than previous theoretical estimates The reason for 

this discrepancy is explained in Fig. ^p. Previous cal- 
culations used a continuum model, which is equivalent 
to including acoustic but not optical phonons. However, 
almost all the binding is due to optical phonons, which 
couple far more effectively to the electrons. If we repeat 
the calculations using only the first three phonon bands, 
we obtain only a weak binding, in excellent agreement 
with continuum calculations 

In conclusion, we find that in semiconducting car- 
bon nanotubes, the scattering of electrons or holes by 
phonons is strongly energy dependent, both at low en- 
ergy and around the optical phonon energies. Neverthe- 
less, the mobility at low carrier densities can be described 
fairly accurately by a simple formula, which should be 
useful in the analysis of nanotube devices. The polaronic 
binding energy is much larger than expected, giving a 
significant renormalization of the bandgap. 
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